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Abstract 

We calculate the one- loop QCD virtual corrections to all helicity amplitudes 
for parton level virtual Compton scattering processes. We include the amplitudes 
both on quark target process 7*g — ► and on gluon target process 7*5 — ► 73. 
The infrared pole structure of the amplitudes is in agreement with the prediction 
of Catani's general formalism for the singularities of one-loop amplitudes, while ex- 
pressions for the finite remainder are given in terms of logarithms and dilogarithms 
that are real in the physical region. 



1 Introduction 



Photons, being real or virtual, are known to be clean probes of the internal structure 
of the nucleon. One of the important processes of photon is Compton scattering, which 
refers to elastic scattering of a photon off a charged object. The Compton reaction 
in different kinematical situations can provide different physical information about the 
behavior of quarks and gluons in the nucleon. There are two complementary kinematical 
regions, which are deeply virtual region and wide-angle region. The region of deeply 
virtual Compton scattering (DVCS) is characterized by small momentum transfer from 
the initial to final nucleon and a large photon virtuality, while in the region of wide-angle 
Compton scattering (WACS) the situation is reversed. 

The theoretical development on DVCS has revealed that in the Bjorken limit, this 
process is dominated by the simple handbag mechanism in which a quark (an antiquark) 
in the initial nucleon absorbs the virtual photon, immediately radiates a real photon, 
and falls back to form the recoiled nucleon QjJ |2l El- It has been proved that DVCS 
amplitude can be factorized into the finite perturbative parts and the nonperturbative 
parts represented by collinearly-divergent terms, which correspond to the matrix elements 
of a class of newly introduced generalized parton distributions (GPDs) jHElEl- These 
GPDs are hybrid objects, which combine properties of form factors and of ordinary parton 
distributions. 

For WACS process, there is a general agreement that in asymptotic large momentum 
transfer region, the amplitude for Compton scattering is dominated by Brodsky-Lepage 
hard scattering picture [7j, which is given by the convolution of a hard-scattering am- 
plitude of collinear constituent partons and the distribution amplitude of hadrons. The 
perturbative contribution has been calculated in 001101 to leading twist accuracy. How- 
ever, the cross sections predicted in that approach are way below the existing Compton 
data, unless strongly asymmetric, i.e. end-point concentrated distribution amplitudes are 
used. This results in the conclusion that the cross sections are dominated by contributions 
from the soft end-point regions, where the assumptions of the leading twist perturbative 
calculation break down. Thus, the hard scattering model of Brodsky-Lepage, although 
likely to be the true asymptotic picture for exclusive reactions, does not seem to be dom- 
inant at moderately large momentum transfer, the kinematical region being accessible to 
present-day experiments. 

Recently, a new mechanism has been introduced, in which the physics of handbag di- 
agram is also of importance even for WACS reaction. As has been argued in [TTlll2[ IT3]. 
at moderately large momentum transfer, Compton scattering off protons approximately 
factorizes into a hard parton-photon subprocess and a soft-proton matrix element de- 
scribed by new form factors specific to Compton scattering. These new form factors 
represent moments of GPDs and can be modelled by overlaps of light-cone wave func- 
tions, which provide the link between exclusive and inclusive reactions. Using light-cone 
wave function overlaps as a model of the GPDs, detailed predictions of cross sections and 
polarization observables for real and virtual Compton scattering have been achieved in 
|12[ ITS] by using the leading order (LO) results for the hard parton-photon subprocess. 
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In addition, the perturbative contributions appeared in the factorization formalisms for 
both DVCS and WACS processes have been calculated to one-loop level so far. The au- 
thors of [5J El El Ej have considered the one-loop corrections to DVCS process, while 
the next-to-leading order (NLO) QCD corrections to the parton level subprocess of real 
WACS processes have been calculated in ^7j. Some numerical predictions were also 
presented in the latter paper J7j based on the GPDs model proposed in [T2"] . 

The purpose of the present paper is to calculate the NLO corrections to the parton 
level virtual Compton scattering in wide-angle region. Our results together with the 
GPDs form factors will provide the necessary ingredients for calculating the wide-angle 
virtual Compton scattering at the NLO. These results may also be used to compare 
with other theoretical results in order to test the handbag mechanism and facilitate the 
interpretation of future experimental data that might be obtained at Jefferson Lab or at 
an ELFE-type accelerator at DESY or CERN. 

The paper is organized as follows: Starting with the kinematical discussion, we cal- 
culate the LO amplitudes in Sect. 2. Detailed calculations about the NLO results are 
given in Sect. 3. Sect. 4 is devoted to the check of our calculations by comparing the 
divergent parts of our formulas with Catani's general formulas, and also by comparing 
the expressions of our formulas in the — t — > (DVCS) and Q 2 — > (real WACS) limits 
with those results appeared in El El El El E] ■ The paper terminates with a few 
concluding remarks (Sect. 5). 



2 LO amplitudes 

For the parton level virtual Compton process 

7*(g) + parton{p) — > 7(9') + parton(p'), (1) 

we work in the center of mass system (CMS) of photon and parton. Neglecting the parton 
mass the different particle momenta can be chosen to be 

<f = (g ,0,0,|p|), 

p" = |p|(l,0,0,-l), 

q'v = |p'|(l,sin0,O,cos0), 

p'A* = |p'|(l, -siiifl, 0,- cos 0), (2) 

in that frame. Here 9 is the scattering angle of the outgoing real photon. q , |p|, |p'| and 
9 are related to the Mandelstam variables s = (p + q) 2 ,t = (p 1 — p) 2 and to the virtuality 
(q 2 = —Q 2 ) of the incoming photon by 



Q 2 , , s + Q 
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(a) 




Figure 1: Feynman graphs for Compton scattering off on-shell quarks, a) is the LO graph, 
the others represent the NLO QCD corrections. Graphs with self-energy corrections to 
external fermions and those with interchanged interaction points of the photons are not 
shown. 

In CMS, the polarization vectors of incoming virtual photon are 

e 3 (±) = T -^(0,l,±z,0), e 9 (0) = ^(|p|,0,0,g ). (4) 

The polarization vector of the outgoing real photon is given by 

£qi{fJ>') = —7=(0,—fjfcos9,—ijfi'siD.9). (5) 

The tree amplitudes only receive contribution from scattering with quark. There are 
two Feynman diagrams which contribute to the LO amplitude for the reaction 7*g — > 7g; 
a half of them is shown in Fig^a). The other half will be taken into account by using 
the crossing symmetry, i.e. the simultaneous replacement of q <-» — q' and e q <->■ e*,. The 
LO amplitude is 

= ey M ^ x , (6) 



where e q is the electric charge of quark q, and /i(A) and //(A') denote the helicity of the 

f (o) 



initial and final state photons (partons), respectively. The amplitude M^\, A is written 
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(b) 



Figure 2: Sample Feynman graphs for photon-gluon scattering. 



as 



M (0) 



u(p', X' 
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Working out in the photon-parton CMS, we explicitly find 
(o) o l~s~s + Q 2 
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-0, 
0. 



u(p,X).(7) 



For the sake of legibility, explicit helicities are labeled only by their signs. Since the 
quarks are taken as massless, there is no quark helicity flip amplitude, i.e. M m '_a )M a = 0, 
to any order of a s . Other helicity amplitudes can be obtained from those given in (JSJ) by 
parity invar iance: 



M_ 



-A',-^-A = (-ir- A ^ +A M M , A , 



(9) 



3 One loop processes 



The NLO corrections to 7*g — > can be calculated from the Feynman graphs (b)- 
(e) depicted in Figd We work in Feynman gauge and use dimensional regularization 
(n = 4 — 2e). There are two mitigating factors which simplify the NLO corrections. 
The first is that the calculation does not contain ultraviolet (UV) singularities since the 
graphs in Fig. do not contribute to the renormalization of the strong, electromagnetic, 
or weak coupling constants. The second is that the self-energy insertions on the external 
quark lines vanish due to the cancellation of the UV and infrared (IR) divergences |19j . 
Basically, what happens is that the UV and IR poles cancel when one does not distinguish 
between them. 
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We evaluate these amplitudes using the Feynman parametrization technique. The 
loop integrals associated with the four-point function from the box diagrams shown in 
Fig^e) are very difficult to evaluate when powers of the loop momenta appear in the 
numerator. However, one can express these tensor integrals in terms of lower rank tensor 
integrals with the same number of propagators and lower rank tensors with fewer prop- 
agators |20j . In the end, the four-point functions with powers of the loop momentum 
in the numerator are reduced to a four-point functions with a constant numerator and 
three- and two-point functions which are easier to evaluate. The all loop integrals can be 
reduced to a set of 13 scalar integrals which are given in Appendix. 

The NLO corrections to the 7*g — > ^yq amplitudes read 



TY 



-ty \2e 2 AeJ A 12 



{2su-tQ 2 )(s + u) + 2(s + t)su (Ast + su + At 2 + 2tu)su Q 2 
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(13) 



a, „ , f 1 / 1 3 \ 9 7T 2 1 / -i 



7T J w ( (-t) £ V2e 2 4e/ 4 12 2 V s 



u s 2 -2ut-u 2 -t f 2 {s,t,u) t \ 
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where the functions /(e) and fi(s, t,u),i = 1, 2, 3 are all defined in the Appendix. We see 
that the NLO amplitudes possess both non-zero imaginary parts and non-zero photon 
helicity flips. 

At the one- loop level, there is a complication which we have to discuss next, namely 
gluons have to be considered as active partons as well. In contrast to the case of quarks, 
the partonic amplitudes now allow parton, i.e. gluon helicity flips to occur. The photon- 
gluon amplitudes can be calculated from the three graphs shown in Fig. EJ There are 
three further graphs contributing to order a s which however reduce to the first three ones 
by reversing the fermion number flow. The fermion loop contributions satisfy the color 
decomposition 

Mf) y ^ x = E^e 2 ^Mj« M , (16) 
i 

where % and a 2 are color indices of incoming and outgoing gluons. After some algebra, 
we find for the gluon helicity non-flip amplitudes 

++ ' ++ 7T \ 2s(t + uy K J (s + u) 2 (t + u) -t 
(2st- su + 2t 2 )u^Q 2 | (s-Q 2 )tu \ 



(s + t) 2 (t + u) -u (s + t)(s + u)(t + u) 
,,„m oc s [ (2t 2 + 2tu + u 2 )s , . . (2t + u)s/ -t 

= i p 2uHt + u) /,(.,«,,) + L_L( ln _ + M , 
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f tQ 2 tuQ 2 Q 2 st \ 
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,f utQ 2 (3s + 2u)tuQ 2 Q 2 

M y + \ _ + = — <^ —r- — -f 2 (s, t, u) - — — ^t— — r ln— 
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(3s + 2t)tuQ 2 _ ln _£^ _ (g - Q 2 )tu \ 
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and for the gluon helicity flip amplitudes 

(23) 

mQ 2 su tuQ 2 Q 2 \ 
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7T 
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a 



uts uQ \/2utsQ 



7T V V 2 t 2 (t + u) J ° K ' ' ' (s + t)(t + u) 

V?ut~sQ ( Q 2 . \ V2uTsQ(st + su + t 2 + 2tu). Q 2 \ 
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(t + u)t \ s J (s + t) 2 (t + u)t -u j 

We find that for all helicity non-flip amplitudes, the scaling behavior of — t — > are 
satisfied, which is demanded by the angular momentum conservation. 



4 Comparison with other results 

We find that those helicity amplitudes which are non-zero at LO, turn out to include 
infrared divergent parts which go as an universal form 



<ui^>(?+!)(^p (29) 



where the 1/e 2 term appears as a consequence of overlapping soft and collinear diver- 
gences. Catani presented a general formula for the structure of infrared divergences of 
any QCD amplitude [21]. For the case of the one-loop amplitude, Catani's formula is 



-e^W ^ 1 / 2 1 l\v-™ m U 2 e~ lX ^\ a s S e , . 

M {0) ; -V^ Ti 2 — + 7i- VT r T, 30 

2Y(l-e)^T 2 \ 1 e 2 h eJf^ t J \ 2 Pi ■ Pj J 2n K } 

where S £ is the typical phase-space volume factor in d = 4 — 2e (je — ~ ^(1) — 0.5772 • ■ • 
is the Euler number) 

S e = exp[e(ln47r + V(l))]- (31) 

In (|3()|1. the sum is over all external partons of the amplitudes. A^- = +1 if i and j are 
both incoming or outgoing partons and Ay = otherwise. The color charge T, = {Xf} 
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is a vector with respect to the color indices a of the emitted gluon, and an SU(N) 
matrix with respect to the color indices of the parton i. Applying (j3*Uj) to the process we 
consider here, one may convert two of the partons to photons by setting T, • T 7 — > and 
Tj • Tj = — T 2 to obtain the simplified formula 

- M (o)gi exp(dn4^) ^ / I + 3 N V ^ 



2tt r(l-e) Ve 2 2e 

the divergent part of which is exactly equal to the one given in eq. ()29j1 . This agreement 
provides a stringent check of our amplitudes. 

To perform further checks on the reliability of amplitudes, we compare the amplitudes 
of real photon-parton scattering processes obtained from our formulas with those already 
existed results. The helicity amplitudes with one-loop corrections to real Compton pro- 
cess 7g — > 7g can be derived from (|10j ) -([15| ) by taking Q 2 = 0, and they are 

q (i) _ t( ^C F a s { 1 /l 3\ 7 7T 2 



+ — In — h -In h — ( In — + 7r ) + inln — > J , (33) 




—u 



Ml {1 2. 



+ 

+ +H r ^ + - ln — ) — . (34) 

= -^h/— + \/— h (35) 

<U = -^IJ— + h (36) 



AfRSV = (37) 

r(i) 
-,o- 



M| ( o } _ = 0. (3* 



We can see that the contributions from longitudinal photon disappear and all other 
amplitudes agree with those given in |17j . Applying Q 2 = to the real photon-gluon 
scattering process jg — > jg, (JIZJ)-(j2BJ) can be simplified as 

't 2 + u 2 (. 2 t 2 \ s + 2tt \ 

— — ln 2 - + vr 2 + ln- + l , (39) 

2s z \ u J s u J 

V + t 2 2 -t t-s-t ft-s s 2 + t 2 , -Al . , 
^_^l n 2_ + i n _ + i + m i + —In— (40 

2u z s us \ u u s \ 



m ++,++ 
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Mill- = ~, (41) 

7T 

M 9 + ( IU = ^, (42) 

7T 

M 9 + % + = 0, (43) 

M£! _ = 0, (44) 

Mf; ++ = ^, (45) 

s t 

In— )| (17) 




S 2 + M 2 , — M 




t 2 

M 9 + {1) _ + = ^, fi> 

7T 

M 9 + % + = 0, (49) 
M£! _ = 0, (50) 



which are also equal to the results obtained in ^7] and 

To compare with the DVCS results, let us take the forward limit, i.e. —t — > 0. Then, 
from (JHJ), we can see that most of the LO helicity amplitudes vanish except 



Since the scalar integrals (|58jl. (|59*|l. (jSOjl and (fTOj) in Appendix are not well-defined at 
the point — t = 0, our methods described above can not be directly applied to the case of 
DVCS. However, this difficulty may be avoided by taking the trace over the numerator 
of box integrals. In the DVCS limit (—t = 0), from (J2J) we have p' = j^qz P and the 
product of two Dirac spinors can be written as u a (p', \)up{p, A) ~ (1 + X-y^p 1 , therefore 
the numerators of the box integrals is of the trace form. The trace of the numerator can 
be expressed as a sum of terms which cancel one of the propagators in the denominator. 
This can be done because k 2 and k ■ p can be written as linear combinations of k 2 and 
(k + p) 2 , where k is the loop momentum, and the trace only includes these two type of 
terms. Since we have shown that the numerator will be linear combination of two different 
denominators, the four-propagator integral will in general become two three-propagator 
integrals. After this procedure, we do not need to carry out the calculation of the four 
questionable scalar integrals mentioned above to obtain the DVCS amplitudes, and the 
calculation becomes much simpler. After some algebras, we find the NLO results for the 
DVCS helicity amplitudes 



M (i) _ M (o) 





In 
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+9-ln^-^(ln^-^-(ln--.vr] }, (52) 

M (1) - M (0) ( { 1 1 (o | 2111—^1 

. \ 2{s + Q 2 ) , w 2 -tt\ 

m J — — V~ ln w (53) 

The 1/e 2 divergent terms vanish and only 1/e divergence is left, this is due to the fact 
that there is no soft divergence in the — t = limit. Except for an overall normalization 
factor, the symmetric and antisymmetric quark amplitudes defined in j3] can be written 

as 

Tf ) = Mj_i_ + M i+J+ , 7f 1 = Mj_ j_ - M j+J -+. (54) 

The CM frame chosen in this paper is different from the CM frame chosen in jS] , and the 
relation between the kinematic variables used in these two papers can be expressed as 

1 ~ XB -Q 2 , Q 2 = 2Q 2 , (55) 



x B 

where Q 2 represents the variable Q 2 defined in After the substitutions of the 
two amplitudes derived from (p>4*j) are exactly equal to the corresponding equations (24) 
and (25) in |5| when the limits of — t —>■ and M 2 —>■ are taken for the later case. By 
checking in the similar way, our results in DVCS limit are also shown to be in agreement 
with other existing one-loop results in [Ti l ITH j ITfij. 



5 Conclusion and discussion 

As a complement to the results given in [T2], we have calculated the NLO QCD corrections 
to the subprocess helicity amplitudes for virtual Compton scattering. The divergent 
parts, being consistent with Catani's general formulas, exhibit the same behavior as 
those in electromagnetic form factor and therefore can be attributed to the form factors 
introduced in [T2] which are specific to Compton scattering. Our finite results can be 
used to the predictions for various virtual Compton observables and are compared to the 
leading contribution given in ^3] • One of the interesting physical quantities is the beam 
asymmetry for ep — > epj 

da{+) - da(-) 
~ da(+) + da(-Y {bb) 

where the labels + and — denote the lepton beam helicity. This quantity depends on 
the relative phase between the complex virtual Compton amplitudes and the real Bethe- 
Heitler ones. In the LO soft physics approach, Al is zero because all amplitudes are real. 
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However, since the a s corrections in the photon-parton subprocess include imaginary 
parts, Ai may become non-zero in the NLO soft physics approach [22] • This must be a 
good test for the handbag mechanism by comparing with the data which will be obtained 
in the forthcoming experiments at Jefferson Lab or ELFE-type accelerators. 

Acknowledgments 

This work started when one of the authors(H.W.H.) was at Kobe university as the post- 
doctoral research fellow (No. P99221) of the Japan Society of the Promotion of Science 
(JSPS) and H.W.H. would like to thank the Monbusho's Grand-in-Aid for the JSPS 
postdoctoral fellow for financial support. The authors also acknowledge useful comments 
from Dieter Miiller. 

References 



[i] 

[2] 

[3] 
[4] 
[5] 
[6] 
[7 
[8 

[9 

[10 

[11 
[12 



F. M. Dittes, D. Miiller, D. Robaschik, B. Geyer and J. Horejsi, Phys. Lett. B 209, 
325 (1988). 

D. Miiller, D. Robaschik, B. Geyer, F. M. Dittes and J. Horejsi, Fortsch. Phys. 42, 



101 (1994) |arXiv:hep-ph/9812448|. 

X. Ji, Phys. Rev. Lett. 78, 610 (1997); Phys. Rev. D55, 7114 (1997) 

A. V. Radyushkin, Phys. Rev. D 56, 5524 (1997) || arXiv:hep-ph/9704207| . 

X. D. Ji and J. Osborne, Phys. Rev. D58, 094018 (1998) |arXiv:hep-ph/9801260| . 

J. C. Collins and A. Freund, Phys. Rev. D59, 074009 (1999) |arXiv:hep-ph/9801262| . 

G.P Lepage and S.J. Brodsky, Phys. Rev. D22, 2157 (1980). 

G. R. Farrar and H. Zhang, Phys. Rev. D41, 3348 (1990) [Erratum-ibid. D42, 3348 
(1990)]. 

A. Kronfeld and B. Nizic, Phys. Rev. D44, 3445 (1991); Erratum Phys. Rev. D46, 
2272 (1992). 

M. Vanderhaeghen, P.A.M. Guichon and J. Van de Wiele, Nucl. Phys. A622, 144c 
(1997). 

A. V. Radyushkin, Phys. Rev. D58, 114008 (1998) larXiv:hep-ph/9803316| . 

M. Diehl, Th. Feldman n, R. Jakob and P. Kroll, Eur. Phys. J. C8, 409 (1999) 
larXiv:hep-ph/9811253| . 



12 



[13] M. Diehl, T. Feldmann, R. Jakob and P. Kroll, Phys. Lett. B460, 204 (1999) 
|arXiv:hep-ph/9903268| . 



[14 
[15 

[16 

[17; 

[18 
[19 
[20 
[21 
[22 



A. V. Belitsky and D. Mliller, Phys. Lett. B 417, 129 (1998) |arXiv:hep-ph/9709379| . 

L. Mankiewicz, G. Piller, E. Stein, M. Vanttinen and T. Weigl, Phys. Lett. B 425, 
186 (1998) |arXiv:hep-ph/9712251| . 

A. V. Belitsky and A. Schafer, Nucl. Phys. B 527, 235 (1998) 
|arXiv:hep-ph/9801252| . 

H. W. Huang, P. Kroll, T. Morii, Eur. Phys. J. C23, 301 (2002) 
jarXiv:hep-ph/0110208| . 

Z. Bern, A. De Freitas, L. Dixon, JHEP 0109,037 (2001) [ arXiv:hep-ph/0109078| . 
W. L. Van Neerven, Nucl. Phys. B268, 453 (1986). 

G. Passarino, M. Veltman, Nucl. Phys. B160, 151 (1979). 

S. Catani, Phys. Lett. B427, 161 (1998) |arXiv:hep-ph/9802439| . 

H. W. Huang and T. Morii, in preparation. 



13 



A Scalar loop integrals 



The loop integrals from the one-loop graphs of Fig. Q and Fig- 121 can be reduced to a 
set of 13 scalar integrals which are given in this Appendix. They all include a common 
overall factor 

»r(l + e)(47r/z 2 ) e 

m = — ^ — • (57) 

Among these 13 integrals, the four-point integrals D l (i = 1 — 3) and the three-point 
integrals Cq (i = 1 — 6) are infrared divergent and ultraviolet finite, while the two- 
point integrals Bq (i — 1 — 4) include only UV divergence. However, the UV divergent 
term containing e cancels when the two-point functions are combined to form the tensor 
integrals and the UV divergences of the individual graphs also cancel in the sum and, 
therefore, the final amplitudes are UV safe. 

The four-point functions that we need can be expressed as 



Dl = ^ 



d n k 1 



(2?r)" (k 2 + ie) [(p + k) 2 + ie][(p + q + k) 2 + ie) [(p' + k) 2 + ie] 



Dl = ft* J 



/(e) f 2 / 1 1 1 \ vr 2 , , J 

d n k 1 



(tt)™ (k 2 + ie) [(p + k) 2 + ie] [(p - g' + A;) 2 + ie] [(p' + k) 2 + ie] 



/(e) f 2 / 1 1 Mi 2 ,, J 

d n k 1 



(2tt)™ (F + ie) [(p + A;) 2 + ie] [(p + g + A;) 2 + ie] [(g' + /c) 2 + ie] 



/(e) f 2 / 1 1 1 \ vr 2 „ , A 

where 

fi(s,t,u) = In 2 +ln — In 21n — In— + Li 2 



s + Q 2 s + Q 2 s s + Q 2 s V s + Ql 

-Ljl + ^-Ljl-^+LiJl + ^+Li ' S 



t y \ sty V sy ys + Q 2 

s + q 2 



-2Mn^— , (61) 



f 3 (s,t,u) = In 2 ^ 9 +ln In 21n^— -In + Li 2 



2 



s + Q 2 s + Q 2 s s + Q 2 s \s + Q 

-Li 2 f 1 + — ) - Li 2 (1 - ^ ] + Li 2 f 1 + + Li 2 ' 



u I \ su ) \ sj \s + Q 



-u 



-2Mn— 2 . (62) 
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2 / —t\ , / — t \ , / S \ I — t 



For Q 2 > -u, 

+Lb fen?) • (63) 

while for Q 2 < —u, 

t t , S t , t 7T 2 1, 



f 2 (s,t,u) = In — — — + ln — — —In 21n — — —In —--In 



u + Q 2 u + Q 2 -u u + Q 2 u 3 2 u + Q 2 

In the real photon case (Q 2 — ► 0), we have the simplified results 

fi(s,t,u) = -In 2 — -2mln— , (65) 



u 

The three-point functions are 

d n k 



h(s,t,u) = -In 2 2i7rln , (66) 

f 2 (s,t,u) = -ln 2 --vr 2 . (67) 



^ J (2tt)« (k 2 + ie)[(p + k) 2 + ie}[(p + q + k) 2 + is] 
/(e) f 1 111 



(s + Q 2 )\(-s-iey (Q 2 ) £ Je 2 ' 
d n k 1 



C 3 = 



(2tt)" (k 2 + ie) [(p> + k) 2 + is] [(p + g + A;) 2 + is] 

/(e) [1 1 

s \ e 2 (— s — ie) e 6 J ' 
<Tfc 1 



(2tt)" (A; 2 + ie) [(p + A;) 2 + ie] [(p' + A;) 2 + ie] 
/(e) [1 1 7T 2 ] 



* \e 2 (-t) e 6 



(2tt)« (A; 2 + ie) [{k - q) 2 + ie] [(jfe - q') 2 + ie] 
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(68) 



(69) 



(70) 



m \ i Lli 

(t + Q*)\(-ty (Q 2 yje 2 > 

d n k 1 
(27r) n (k 2 + ie) [{p + k) 2 + ie] [{p - q> + k) 2 + ie] 

/(c) fl 1 



o — /i 



7T 

■u [ e 2 (— w) e 6 
d n k 



C e = „* [ dnk \ 

^ J (2tt)« (k 2 + ie){(p> + k) 2 + ie][(p' - q + k) 2 + ie] 



m \ 1 Lli 

(u + Q 2 )\(- u y (Q 2 yje 2 - 



The two-point functions are 
Bl = 



cTA; 1 
(27r)« (A; 2 + ie) [(p + g + k) 2 + ie] 



5o 2 



R 3 — ,, 2e 



'1 1 



(71) 



(72) 



(73) 



^7(^ + 2 J' < 74 > 

^ 7 (2tt)™ (A; 2 + ie) [(p - p> + k) 2 + ie] 

M (^ + 2), (75) 

eTfc 1 
(27r) n (fc 2 + ie)[(g + A;) 2 + ie] 

= m ^m +2 )' (76) 

cf 1 *; 1 
(2tt)™ (k 2 + ie)[{p - q' + k) 2 + ie] 



e 
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